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Abstract 

We investigate the two main uncertainties on the extraction of the nonper- 
turbative matrix elements A and Ai from moments of the inclusive rare FCNC 
decay B Xgj. The first one is due to unknown matrix elements of higher 
dimensional operators which are estimated by varying the matrix elements 
in a range as suggested by dimensional analysis. The effect of these terms 
is found to be small. A second uncertainty arises from a cut on the photon 
energy and we give model independent bounds on these uncertainties as well 
as their estimates from a simplified version of the ACCMM model. 
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I. INTRODUCTION 



Heavy Quark Effective Theory (HQET) has proven to be a useful tool to calculate 
inclusive decays of hadrons containing one heavy quark P,^. The presence of two widely 
separated scales, the mass of the heavy quark, mq, and Aqcd, allows us to perform an opera- 
tor product expansion (OPE) on the inclusive decay rates. At leading order this reproduces 
the parton model result with higher order terms suppressed by powers of Aqcd/'^q- To 
order (Aqcd/'^q)^ they are parametrized by 

Xl = ^{Hiv)\Ki^D)'K\Hiv)) 

= ^^(H{v)\K^{-za,,)G^''K\H{v)), (1.1) 

and 

Here is the spinor of a heavy quark with velocity v and pole mass mg, M is the mass of 
the heavy meson and dn = 3(— 1) for pseudoscalar (vector) mesons. The gluon field strength 
is defined by Gf^u = [iD^, iD^] and the covariant derivative by = d^ — igT'^A'^^. Knowledge 
of the values of these parameters is important for a reliable extraction of the CKM matrix 
elements Vcb and Vuh- The matrix element A2 can be extracted from the observed mass 
splitting between the vector and the pseudo-scalar mesons, 

A2= ^^^^ ^0.12GeV^ (1.3) 

whereas no such simple relation exists for Ai and A. 

Many attempts have been made to extract values of these parameters. Two model 
independent methods have been suggested using moments of the electron energy spectrum 
^ and the invariant hadronic mass spectrum in the inclusive decay B XcCv. In both 
approaches the moments are proportional to linear combinations of the two unknown matrix 
elements. Unfortunately, for the electron energy spectrum, the two moments are given by 
almost the same linear combination of A and Ai, thus, even small corrections to the moments 
lead to large uncertainties on the parameters 0. For the invariant mass spectrum, the two 
moments determine two almost orthogonal linear combinations, but corrections suppressed 
by {^Qco/^bY huge for the second moment. This yields uncertainties for A and Ai 
comparable to those from the electron energy spectrum P] . 

In Ref. [0 it has been shown that the mean photon energy and its variance in the inclusive 
decay B — ^ ^^7, which are given to order 1/m^ by 

_A^A AW 

^ ^' 2 \ 2ml I 2 ' 

var(E,) = (1.4) 
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where var(£'^) = {E"^) — (£"7)^, determine the two parameters independently, thus, this anal- 
ysis yields two orthogonal bands in the A-Ai plane. In this paper we study the theoretical 
uncertainties on the extraction of Ai and A from this decay. Besides perturbative corrections, 
which we do not consider in this paper, there are higher order non-perturbative corrections 
and we estimate their effects by calculating the l/m\ corrections to the first two moments. 
The unknown matrix elements of the dimension seven operators are varied in their expected 
range ±Aqqj-,. Constraints on some of the matrix elements can be obtained from the known 
mass splitting of vector and pseudoscalar mesons and from a vacuum saturation approxi- 
mation. They will be discussed in the following sections. In addition, the measurement of 
the differential decay spectrum involves imposing a cut on the photon energy close to its 
endpoint |jlO|. As the OPE breaks down close to the endpoint of the photon spectrum, we 



anticipate errors which depend on this cut. We will give model independent bounds on those 
uncertainties as well as estimates using the ACCMM model |]lT] . 



II. THE OPERATOR PRODUCT EXPANSION AND MATRIX ELEMENTS 

The differential decay rate of the inclusive decay B — > is given by 

= ^^d[p.s],um^b){i\R\m{q), (2.1) 

where 

W{q) = J2 d[P.S.]hadr{27rYS\PB ~Px- q){B\0^\Xs){X,\0\B). (2.2) 

R and O are the radiative and hadronic part of the operator that mediates the decay, 
respectively. The matrix element of R can be calculated perturbatively whereas W{q) is 
related to the forward scattering amplitude 

T(g) = {B\T{0^0}\B) (2.3) 

via the optical theorem 

W{q) = 2lmT{q). (2.4) 

As has been argued in @J^, the time-ordered product can be expanded by performing an 
OPE 



1 

rrib 



1 1 

h 



(2.5) 



The first three terms in this expansion have been calculated for an arbitrary hadronic oper- 
ator of the form O = bTs in Ref. [Q. We have extended this result by calculating the 1/mf 
contribution to the time-ordered product, O3. The general expressions for these operators 
are quite lengthy and we present them in Appendix ^. 

To obtain T(g), we have to determine matrix elements of these operators. Following Ref. 
||r^ closely, the matrix elements at leading order. 
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{B\bTb\B), (2.6) 

are only nonzero for F = 1 or F = 7^^. The matrix element of the conserved vector current 
is given by 

{B\h'b\B) = 2P^, (2.7) 

where = Mbv^ is the momentum of the heavy B meson. Because of this relation the 
dimension three operators are most conveniently written in terms of the full QCD spinors 6, 
rather than the effective spinors h^. At higher orders in the 1/m;, expansion the spinors are 
expanded and written in terms of h^. Because of Lorentz invariance, the matrix elements of 
the dimension 4 operators can be written as 

{B{v)\KTD^h,\B{v)) = v^{B{v)\KT{iv ■ D)h,\B{v)), (2.8) 

and vanish at leading order due to the equation of motion of the heavy quark field 

{iv ■ D)K = 0. (2.9) 

At higher orders in the l/rrib expansion corrections to this arise which can most easily be 
incorporated by using the equation of motion to order l/rrih 

Kiiv ■ D)K = -^K{iDfK + -^Kia^.G^'^K + O (W . (2.10) 
Zmi) 4mb / 

Matrix elements of general dimension five operators can be parameterized by 

{B{v)\KViD^iD,K\B{v)) = MsTr |fP+ {^\i{g^. - v^v,) + ^As^a^.) P+} , (2.11) 

where P+ = ^(1 + |^) is the projector which projects onto the effective spinor h^. Finally, the 
dimension six operators can be parameterized by the matrix elements of two local operators 

-^^{B{v)\K{iDa){iD^){iDi3)h^\B{v)) = \pi {g^p - VaVf^) v^, 



2Mb' ' " ' " " ' ' " 3' 

{B{v)\K{iDa){iD^){iDp)-islbhv\B{v)) = -p2ievaf3sv''v^ (2-12) 



2Mb' ' " ' " 2 
and by matrix elements of two time-ordered products 

-{B{v)\K{iDfh,t / d^x / Ci{x)\B{v)) + h.c. ' ' 



2Mb J J-00 nib 

-—i^B[v)\K-{-ia^,)G^-Ki \ d'x / Cj{x)\Biv)) + h.c. = \ (2.13) 

2Mb 2 J J-00 rrib 
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III. HIGHER ORDER CORRECTIONS TO THE DIFFERENTIAL DECAY RATE 



Integrating out the W boson and the top quark the FCNC transition b s is 
governed by the effective Hamiltonian 



n 



eff 



4Gf 

"71" 



(3.1) 



In the leading logarithmic expansion, the decay B Xg'-f occurs only through the operator 



O7 



1677 ' 



(3.2) 



which leads at the parton level to the total decay rate 

^ Gj,\Vt,V:faCM' 



3277^ 



ml. 



(3.3) 



Higher order terms in the OPE arise from the expansion of the s-quark propagator in Fig.|l| 



mhf + jt — ^ — rus mfe A 

where A = 1 — 2f ■ g — m^, as well as from the expansion of the full spinors b 

/ ,11. p \ 



(3.4) 



1 + 



2mb 8m? 



+ 



(3.5) 



Starting at order l/mf, additional contributions arise from the difference of the B meson 
states in the full and the effective theory [0. At order 1/ml they are parameterized by 



2M, 



B 



B(v) 



hJiDYK 



6M 



B 



B(v) 



B{v. 
Biv) 



Ai + 



A2 + 



Ti + 3r3 
% + 3r4 

3mfe 



(3.6) 



where the matrix elements Tx-Ts, have been defined in ( |2.13| ). 

To calculate T(g) for the operator Oy, we now use the expressions for the generic oper- 
ators given in Appendix 0, with Fi = PiCTtMu and F2 = o-p^Pn contracted with q^q" , the q 
dependence of the matrix element of R. Parameterizing the matrix elements as outlined in 
the previous section and making the replacements 



Ai ^ Ai + 



A2 + 



nib 
% + S% 
Srrib 



(3.7) 



to incorporate (|3.6|) we find 
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T{q) = -16MBmb{v ■ qf 



1 

A 



Ai /5 — 



2ml 



3A3 



+ 



A, 



2ml 



' 2fd{v ■ qf - 2Av q + 9\ ps (2vq + l 

2m\ 
?,% 1 



2mf 



leading to the total decay rate 
^ ^ ( Ai - 9A2 




A3 



A2 



(3.8) 



1 Ipi - 27p2 Ti + 3T2 - 3 (T3 + 3T4) 



6mf 



+ 



2ml 



(3.9) 



Here Fq is the parton level decay rate as defined in (|3.3|) . To order 1/ml this is in agreement 



with the expressions obtained in Ref. fl^JlSf . 

From the expression of T[q) ( p.8|) we can also obtain the differential decay rate 



dV 

dx 



V [5(1 -x) + A^5\l -x)+ A25"{1 -x)+ A^5"\l - x)] 



(3.10) 



where 



^3 



Ai + 3A2 5pi - 21p2 Ti + 3T2 + T3 + Z% 



2ml 



6ml 



2ml 



Ai 2pi - 3p2 Ti + 37^ 



6ml 



6ml 



6m? 



Pi 



18mf ' 



(3.11) 



and X = 2E^/mb- Again, to order 1/ml, this is in agreement with the expression found 
in Ref. |]12[. We note that none of the coefficients of the dimension seven operators is 
anomalously large compared to the ones of the dimension six operators. This indicates that 
the 1/ml corrections give rise to only small uncertainties on the moments of the photon 
energy spectrum. This will be investigated in much more detail in the next section. 

Recently, Voloshin has shown the existence of a correction to the total decay rate arising 
from a local operator suppressed by powers of AQQj^/ml, which is given by |l^ 



AF = -Fr 



A2C'2 



QmlC-r' 

Cj and C2 are the coefficient of the operators O7 defined in (|3.2| ) and 



O 



2 = SL^J^Lf.CLf.'Jf.CL^, 



(3.12) 



(3.13) 



respectively. At the scale fi = mb the values of these coefficients are C2 = 1.11, C7 = —0.32 
and the correction to the total rate is only about 2.5%, of the same order as the corrections 
due to the leading 1 / corrections. Since the leading terms in the expressions for the mean 
photon energy and its variance are already suppressed by powers of 1/ m^, this operator could 
have a significant effect on their values. Fortunately, the terms Ai-A^, which contribute 



6 



to the higher moments of the spectrum, arise from multiple poles (1/A)" and are therefore 
suppressed by additional powers of 1/mfe. The leading contribution due to this new operator 
are suppressed at least by l/(mbm^) and up to this order we find 



Ar = ro 



A2C2 PlC*2 I3P2C2 



P2C2 



AA2 = AA3 = 0. (3.14) 

These terms are of the same order as the contributions from dimension seven operators 
( p.ll| ) and will be included in the determination of the uncertainties on the extraction of 



the parameters A and Ai. 



IV. UNCERTAINTIES ON A AND Ai 

To order the mean photon energy and its variance can be obtained from ( p.lOD , 

( CT ) and (|3ll ) 

, ^ , Mb- A 13pi - 33p2 H + 3T2 + T3 + 3T4 P2C2 
(^7/ = n 



24ml 8"^^ ISm^Cy 



/c. N -^1 2pi-3p2 T1 + 3T2 

var Lfc^ = . 4.1 

^ 12 12mfe 12mfe ^ ^ 

The size of the matrix elements pi, p2 and Ti-T^ are unknown, but by dimensional analysis 
we expect them to be of order iAq^j^. Based on a vacuum saturation approximation pi is 
assumed to be positive. 

The relation (|1.2| ) for the mass splitting between the meson mass M and the quark mass 
mg also has to be extended to include the l/mg corrections. To this order it is given by 



, T \i + dH\2 , pi + dHp2 Ti + + du {T2 + 

M = mQ + A \ — 2 z — 2 • (4-2) 

2mQ ArriQ 4m^ 



Using the AMb = Mb* — Mb and the AMd = M^* — M^ mass splitting, the relation 

^ ^ _ kK)M| AMb(Mb + A) - MIAMd{Mb + A) 
~ ' ~ ' " M^ + A-«:K)(M^ + A) ' ^^-^^ 

where K,{mc) = (a(mc)/a(mfe))^/^", can be used to eliminate one of the 6 unknown parame- 
ters. 

In order to compare the uncertainties on the values of A and Ai from dimension seven 
operators in this analysis with the corresponding ones in the semileptonic decay B — »• XcCu, 
we estimate the effect of these higher dimensional terms in the same way as in 0J^. Since 
the current measurement of the photon energy spectrum does not constrain the first two 
moments to any reasonable accuracy, we use the hypothetical data (E^) = 2.45 GeV and 
var(£'^) = 0.13 GeV^ in the following analysis (these are obtained from the central values of 
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Ai and A, extracted from inclusive semileptonic B decays, by using Eq. ( p..4|) ). With this 
input, we now use Eq. (|4.1|) to extract A and Ai and determine how their values are affected 
by the dimension seven operators by randomly varying the magnitude of the parameters 
Pi, p2 and T1-T4 in between and (0.5 GeV)^ and imposing the positivity of pi as well 
as the constraint ( [4.3| ). In Fig.^ we show an ellipse in the A - Ai plane which is centered 



about the mean of the distribution and which contains 68% of the points. This should give 
a reasonable estimate of the theoretical uncertainty in the extraction of A and Ai due to 
higher order corrections. As a comparison we also show the corresponding ellipse obtained 
from the hadron invariant mass spectrum of the inclusive decay B XcCu [[7|. Of course, 
the relative position of the two ellipses has no meaning since we have not used experimental 
values of the moments for the photon spectrum. Only the relative size can be compared. 
The standard deviation of the two parameters is given by 

AA^O.OSGeV (4.4) 

and 

AAi^O.OSGeV^ (4.5) 

These uncertainties are much smaller than the expected size of the matrix elements, thus, 
the corrections should be well under control. 



V. UNCERTAINTIES FROM THE PHOTON ENERGY CUT 

A cut on the photon energy to suppress the strong background from other B decays 
reduces the amount of phase space considerably. In its measurement of the inclusive photon 
energy spectrum |TD|, CLEO imposed a cut at 2.2 GeV, about 500 MeV below the maximum 



photon energy. The OPE is only valid if smeared over a region much larger than Aqcd, but 
it has been shown that for a smearing region of order Aqcd the most singular terms 



in the OPE, which are formally of the same order, can be resummed in the so called shape 
function. For an even smaller smearing region of order Aq^j-j/mf, all the subleading terms 
become equally important. 

Since the OPE gets worse the closer the photon energy cut gets to the endpoint of the 
spectrum, one expects uncertainties that depend on the value of this cut. The expression 
for the differential decay rate obtained by performing an OPE contains only delta functions 
and their derivatives which contribute at = mi,/2. A photon energy cut affecting the 
lower bound on the range of integration does not change the results for the moments as long 
it is below = mh/2. The effect of the cut is hidden in the fact that the OPE for the 
differential decay rate gets worse as the smearing region diminishes. The experimental value 
of the higher moments of the photon energy spectrum measured in the presence of a cut 
will be larger than the value which determines the parameters A and Ai as defined in ( p.. 41 ). 
Therefore, a cut will shift the values of the two parameters extracted in this way, leading to 
uncertainties on their values. 

In this section we will take two different approaches to investigate the size of those 
uncertainties. The first will use the shape function analysis to yield model independent 
bounds on the uncertainties and in the second we will estimate them using a simple version 



of the ACCMM model 11 . 
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A. Model independent bounds on the uncertainties 



The nth moment {n > 1) with a cut on the photon energy Eq is given by 



J,^--e{E-Eo)E-§dE 

where 

r^o = / ^dE (5.2) 

Jeo dE ^ ' 

is the total decay rate with a cut. The positive powers of E in the integrand of the numerator 
of Eq. (|5.1| ) weight the higher energy part of the spectrum more, therefore 

M„^o > (5.3) 

where M„ is the moment without a photon energy cut. 

In [|18| a method has been proposed to obtain model independent bounds on the total 
decay rate. In this approach one replaces the step function with a smooth function P{E, Eq), 
obeying 

P{E, Eo) < 1 for Eo<E< E^ax 

P{E, Eo) <0 for < ^ < Eq, (5.4) 

to find the inequality 

YEo>Yp= PiE,Eo)^dE. (5.5) 

JO dE 

The resulting bounds on the total decay rate are shown in Fig.^. One can see that depending 
on the values of Ai and rub, a cut at 2.2 GeV can have a significant effect on the measurement 
of the total decay rate. This is in agreement with the analysis done using the ACCMM model 
19| which has been used by the CLEO collaboration |T^. For a monotonically increasing 



function P{E,Eq), this analysis can also be used to obtain a bound on the moments as 
measured by experiment 

p(E,Eo)E-^%dE p p , , 

J- p 

The most singular terms can be resummed into the shape function 

where the coefficients An now only contain the leading term in the 1/mb expansion. Using 
this we find 
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^\E=mt/2 


(mi, 


yAr,{P{E,Eo)Y% 


E=m,,/2 



If P is a polynomial of order k, only the first k terms in the shape function ( ^.71) contribute. 
In the following analysis we will use the polynomial 

P(E,Eo) = l-f|^^V Eo<E,. (5.9) 



El — Eq 

In order to have a monotonic function which satisfies the requirements ( |5.4| ), Eq has to 
coincide with the value of the photon energy cut and Ei with the maximum photon energy, 
E 

Since we are working to order we have to use a polynomial of order (3 — n) to 

obtain a bound on the nth moment. For the first moment we find using Eqns. (|5.8| ) and 
T9D 



j^Eq ^ 3 (2 - nib) nih (2 - 4 + mf,) + (-4 E^ + S nib) \i + pi 
' - 2 {3 i2Eo-mb) i2Eo-AE^ + mb) + Xi) 



(5.10) 



From this result it is trivial to obtain a bound on the uncertainty on A. A plot of this bound 
for pi = (0.3 GeV)'^ and two different sets of values Ai and nib is presented in Fig.||. For a 
cut at 2.2 GeV the bound on the uncertainty on A is between 100 MeV and 400 MeV, thus, 
it can be of the same order as Aqcd, leading to large errors. In order to be able to determine 
a precise value of A, the cut would have to be lowered to around 2 GeV. Here, the bound 
on the uncertainty is between 20 MeV and 100 MeV. Since Ai is related to the variance of 
the spectrum, which is the difference of two moments, the P-function analysis can not be 
used to bound its uncertainty. This is due to the fact that the variance is sensitive to the 
tails of the spectrum. 

An inequality such as derived in ( |5.5| ) holds true only for a smooth, positive definite 
shape function. In general, there is nothing that guarantees the validity of this assumption, 
although models such as the ACCMM model |TI[] predict a positive shape function. It should 
also be pointed out that the inequality ( p.5| ) is true rigorously only for smooth functions 
dVs/dE. In this approach we use a singular expansion of the differential decay rate, so 
the convergence of this expansion is essential. As said above, the shape function analysis 
is only valid if the smearing region is of order Aqcd, a scale which also sets the size of 
matrix elements of higher dimensional operators. We therefore expect the bounds to break 
down for values of the matrix elements which exceed the available phase space considerably. 
The breakdown of the inequality (|53| ) can be seen by using Em = 2.6 GeV, Eq = 2 GeV, 
Ai = —0.2 GeV^ and nib = 4.8 GeV. Inserting those values into (|5.5|) we find 



0.193-^ + 0.917= ^ < 1 (5.11) 
GeV^ r - ^ ^ 

which breaks down for pi > (0.76 GeV)^. 

Since we could not obtain a bound on the uncertainty on Ai and we are in no position to 
rigorously justify the assumptions we made in this chapter we will now use a simple model 
to calculate the effect of the photon energy cut on the first two moments. 
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B. Uncertainties using the ACCMM model 



We will use a simplified version of the ACCMM model to estimate the value of the 
moments as measured in the presence of a cut on the photon energy. Assuming a Gaussian 
distribution of the relative momentum of the b-quark inside the B meson, neglecting the 
mass of the s-quark and the momentum dependence of the b-quark mass, the spectral 
function is given by [0 

if . 1 exp|-(^|, (5.12) 

where g\ = — y^. 

The difference between the first moment in the ACCMM model with and without a cut 
gives us an estimate of the uncertainty on the value of A and the result is shown in Fig.|^ 
for different values of ge- For a cut at 2.2 GeV, the uncertainty on the parameter A in this 
model is between 20 MeV and 180 MeV, depending on the width of the spectrum. Since 
the real effect of the photon energy cut could easily exceed this estimate by a factor of two 
or three, a cut on the photon energy at 2.2 GeV could destroy the possibility of accurately 
determining the value of A. If the cut could be lowered to ~ 2 GeV, then an accurate 
extraction should be possible. This is in agreement with the model independent results 
obtained in the last section. 

The result for a similar calculation for the variance of the spectrum is shown in Fig.^. For 
a cut of 2.2 GeV, the uncertainty on Ai in this model is between 0.05 GeV^ and 0.3 GeV^, 
Again, considering the fact that this is only a model calculation and it might underesti- 
mate the effect considerably, this indicates that an extraction of Ai from the present CLEO 
measurement might be unreliable. Lowering the cut to ~ 1.9 GeV should enable a precise 
determination of this parameter from the decay B Xg'j. 



VI. CONCLUSIONS 

We have analyzed the uncertainties on the extraction of the two nonperturbative matrix 
elements A and Ai from the mean photon energy and its variance in the inclusive decay 
B Xsj. Besides perturbative corrections which we have not considered here, uncertain- 
ties arise from matrix elements of higher dimensional operators which are suppressed by 
additional powers of Aqc^ /rrih. We have calculated the first two moments up to order l/mf 
and estimated the effect of the unknown matrix elements by varying their values in their 
expected range of magnitude. The resulting uncertainties on A and Ai are given by 

AA^ 0.03 GeV (6.1) 

and 

AAi^O.OSGeV^ (6.2) 

A more serious uncertainty arises from the effect of a cut on the photon energy which has 
to be imposed in order to reduce the large background from other processes. The differential 
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decay rate as calculated via the OPE is given in terms of a delta function and its derivatives 
contributing at = mi,/ 2. The effect of a cut on the photon energy therefore does not 
affect the results for the moments. It is hidden in the fact that the OPE breaks down as the 
cut approaches the endpoint of the photon energy. We have used an approach suggested in 
||T8| to obtain a model independent bound on the uncertainty on A, whereas no such bound 
could be derived for the uncertainty on Ai. The bound indicates that an accurate extraction 
of A is definitely possible for a cut at 2 GeV, whereas for the present cut at 2.2 GeV the 
errors might be large. We have also used a simplified version of the ACCMM model to 
estimate the effect the photon energy cut. The uncertainties depend strongly on the width 
of the spectrum in our model. We again find that an accurate determination of A should be 
possible if the cut can be lowered to about 2 GeV, Depending on the width of the spectrum, 
the cut has to be lowered even further to allow a precise determination of Ai.. 
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APPENDIX A: THE GENERAL EXPRESSIONS FOR THE FIRST THREE 

OPERATORS IN THE OPE 

In this appendix we will present the first three terms in the OPE for general operators 



T{bTis,sT2b} = — 
nib 



Zruh 



Am 



imt 



(AI) 



With the conventions 



D, = d,-igT'^Al 



(A2) 



they are given by 

Co = --bT,{f,-i + ms)T2b 

X 

01 = --hT^-f''T2^D^h + ^{v - qThT^{i) - ^ + m,)T2iD^h 

16 - " 

02 = --^{v - qr{v - qfhV^ii, - i + ms)T2tD^tD^h 



X 

+^hTi{ij-i + ms)T2{iDfh 



x^ 



+ - 



x^ 



qfhTi-i''T2{iDJDp + iDpiD^)h 



-—mshTiia^0T2G''^h 



(A3) 
(A4) 



(A5) 



x" 



12 



X 

--^(7'^ri7"r2 + T^-fl'T^rYDpiD^h 

X 

+^{v - qThV^ii) -i + ms)V2i^iDJDph (A6) 

X 

32 

'v - qY{v - qYhVi-i'^T2{iDJDpiD^ + iDpiD^iDa + iDpiDJD„)h 



3a;3 
32 



3x3 
64 

4 



3x4 
16 



3x3 
16 



V - qyliT^^i} -i + rhs)T2{iD^{iDf + [iDfiD^ + iDpiDJD^)h 

V - qY{v - qY(v - qYhVriij - j + ms)T2(iDJDpiD^ + iDpiD^iD^ + iDpiDJD^)h 

V - qY{v - qYhT,-f^V2{iDpG^p)h 
- q)fhV,{ij - i + ms)T2{iD^G,p)h 



3x3 
16 - 

+\rhs{v - qYhV^i(T'^PT2{iD^, G^p}h 

X 

+^^£/3p-(^; _ q)^^^ _ g),/iri7,75r2{^Z^\ G^p}h 

X 

+-v''v^hV^^-V2{iD^G^p)h 

X 

--v'^v^hT^-i''T2iaP''{iDp, G,p]h 

X 

--^{i-v qY'^hrs - i + m,)r2(iD^G^^)h 

X 

+^{l-v q)v''hV^{i} -i + rhs)T2iaP''{iDp, G,^}h 
--v'^hViY^2l''iDpiDJD^h 

X 

--hTi^''T2iDJiDfh 

X 

3- 

+-hri-f''r2iaP''iDJDJD„h 

X 

+^v''{v - qfhViii^ - ^ + m,)r2YiDpiDJDph 

X 

+^{v - qYhVi^f -i + ms)T2iDa{iDYh 

X 
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-{v- qThV^ij) -i + rns)T2i(Tf^iDJDpiD^h 

--v"hYTijT2iDpiDaiD^h 

X 

X 

+-hia"^V^Y^2iDJDpiDph 

X 

+ — {v - qYhV^ii) + ms)V2{iD)HD^h 

X 
X 

+\{v - qThYViiij -i + rhs)V2-i^iDpiDJDph 

X 

+\{v - qfJiYTi-t'^V^iDpiiD^iDp + iDfszD^)h 

X 

X 
1 a. 

- — {v - qY{v - qYTiYT^if - i + ms)T2iDpiD^iDph 

X 
X 
X 

+\{v - qfJiV^Y^^YiiDaiDp + iDpiD^)iDph 

X 

+^hTi{^ -i + ms)T2Y{iDfiDph 

16 - - 

- — {v - qT{v - qfhV^if - i + m,)r2YiDaiDpiDph 

X 

-^mJiV^ia'^PV^YGapiDph 

X 
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FIGURES 




I S 




I 

FIG. 1. The Feynman diagram for the forward scattering amplitude. 




0.3 0.35 0.4 0.45 0.5 0.55 0.6 

A (GeV) 

FIG. 2. The estimated size of the uncertainties due to l/m^ corrections on Ai and A from 
B Xg'j (solid) and from semileptonic decay (dashed) using the method described in the text . 
The position of the two ellipses has no meaning, only the relative sizes. 
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FIG. 3. The model independent bound on the total decay rate. 
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FIG. 4. The model independent bound on the uncertainty of A. 
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FIG. 5. The uncertainty on A due to a photon energy cut from the ACCMM model. 
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FIG. 6. The uncertainty on Ai due to a photon energy cut from the ACCMM model. 
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